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We show that the Pu{oj) oc uj^"^/^ shear velocity power spectrum gives rise to a 
Pq{uj) oc uj~^/^ power spectrum for a passively advected scalar, as measured in experiment 
[K. Sreenivasan, Proc. R. Soc. London A 434, 165 (1991)]. Applying our argument to 
high Rayleigh number Rayleigh Benard flow, we can account for the measured scaling 
exponents equally well as the Bolgiano Obukhov theory (B059). Yet, of the two explana- 
tions, only the shear approach might be able to explain why no classical scaling range is 
seen in between the shear (or B059) range and the viscous subrange of the experimental 
temperature spectrum [I. Procaccia et al, Phys. Rev. A 44, 8091 (1991)]. 
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1 Velocity shear spectra 

The classical scaling laws in fully developed turbulence are the Kolmogorov [|^] and the 
Obukhov-Corrsin [Q] scaling relations 

E^{k) ~ ej/'k''/', Ee{k) ~ 606.~V3^-5/3 

for the velocity spectrum Euik) and the temperature spectrum Ee{k) of a passively ad- 
vected temperature field. Here, e„ and ee are the energy and thermal dissipation rates, 
respectively. Eqs. (|I|) follow from dimensional analysis. They can also be derived from the 
assumption that the energy flux T„(A;) ~ kv?{k) of the conserved total energy / v?'{x)dx/2 
and the corresponding flux TQ{k) ~ ku{k)Qi^{k) of the conserved total temperature in- 
tensity / Q'^{x)dx do not depend on fc, i.e., T^ik) ~ and T^i^k) ~ £9. Here we have 
assumed that both velocity and temperature fields are stirred only on the largest scales 
Q and, of course, for these large scales we have neglected viscosity and diffusivity. u{k) 
and 6(fc) are the discrete Fourier transforms of the velocity and the temperature fields 
(having dimensions of length/time and temperature, respectively). 

In anisotropic flow, a correction to (0) is expected. Lumley was the first to suggest 
the spectrum 

E^{k)r^eJ/'k-'/'il + aJ^] ' ) (2) 




for an anisotropic flow, e.g., a shear flow. Here, a„ is an unknown dimensionless parameter 
and ks the crossover wave vector separating the shear dominated range (SDR) and the 
classical scaling range. Eq. (|^) has also been obtained by dimensional analysis in terms 
of Clebsch variables [§, 0]. 

In what follows we will be interested in situations where the second term in (|^) domi- 
nates. This definitely is the case in cross spectra Ei2{k), which probe anisotropy. Indeed, 
Ei2{k) oc was measured in boundary layer flow 0. Also, Maloy and Goldburg's 
scaling of the velocity structure function D^'^\r) oc r^/^ in Taylor-Couette flow with an 
oscillatory inner cylinder might possibly be understood in terms of (^ with a large 

We also mention that Eu{k) oc scaling is expected in helical flow, i.e., flow for 

which the helicity / u ■ V x udx is non zero. The corresponding flux is then of order 
Tneiik) ~ k'^v?{k). Assuming that Tueiik) = e^ei does not depend on A; in a hehcity 
dominated range, we obtain 

Eu{k) ~ eHei"'k-'l\ (3) 
The crossover kg to classical scaling (|l]) can be estimated as kg = enei/^u- 



2 Temperature shear spectra 

How does the spectrum of a passive scalar (e.g., the temperature) look, if advected by a 
shear velocity field? From the conserved flux TQ{k) ~ ku{k)Qi^{k) ~ ee and the shear 

different stirring force is discussed in |^. 
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velocity scaling u{k) ~ e„^/^fcy^fc we readily conclude Q{k) ~ e0-^/^e„ ^^^k^ ^/^k -^/^ 
or 

Ee{k) ^ eeeu-'"k^"'k~"\ (4) 
Considering also the high-Zc range of the spectrum as in we may write 

Ee{k)^eeeu-"'k-^"^ + ae[^^^ '^'j . (5) 

The dimensionless constant cte will be estimated below. 

Experimental support for our suggested scaling law @) is given by Sreenivasan [|T0[| , 
who measured the power spectrum Pq{uj) for temperature fluctuations in the wake of a 
heated cylinder for Re\ = 200. The power spectrum Pe{uj) corresponds to E^ik) via 
Taylor's hypothesis [[T^. Indeed, P@{uj) oc uj~^/^ for more than 1.5 decades, see fig. 1. 
The long scaling range is surprising for the small Re\ as it was in . Unfortunately, the 
velocity spectrum in this experiment was not measured. For the cross spectrum one 
would expect Ei2{k) oc as in the boundary layer [§. We thus felt justified to assume 
u{k) oc here, as the flow behind the cylinder is strongly anisotropic. 

Where should we expect the crossover between the SDR and the K41 range? We 
estimate the crossover frequency Us directly by the shear s. With the data given in [p!o[| , 
we find ujs = s = U/L = 885sec~^. Here, U = 1770cm/s is the flow velocity and L = 2cm 
the diameter of the heated cylinder. Our result log]^Q(to'sSec) = 2.9 is indeed very near 
the end of the Pe{uj) oc uj~^^^ scaling range, cf. fig. 1. The K41 subrange in between the 
SDR and the viscous subrange (VSR) seems to be small, if it exists at all for this small 
Rex. Yet for larger Rex it might show up in between SDR and VSR. Along this line of 
argument we can understand why the (effective) scaling exponent of Pe(co') approaches 



the classical value —5/3, when Rex is increased, see fig. 6 of ref. |]T0 . 

To be more quantitative, we compare the Rex dependence of kg and that of the viscous 
dissipation cutoff k^. The latter is calculated to be |jl2| k^^ ^ ISLRex'^^"^ ■ For the former 



we estimate k~^ < L, as in this geometry shear is induced in the flow at some large 
length scale. So for increasing Rex, a classical scaling range can develop between the 
SDR k < a^kg and the VSR k > kd, which dominates the spectrum at large k. From the 
experimental observation that the K41 subrange is still missing for Rex = 200 |10|, we can 
estimate to be = {k(i{Rex = 200)/A;s)^/^ ~ 5. The classical scaling range will grow 
as kd/{ci%ks) = {Rex/^OOY^'^. For Rex = 1000 it is already larger than one decade and 
more or less classical scaling exponents will be extracted from the temperature spectra, 
which accounts for the measurements summarized in fig. 6 of p!0|. 



3 Rayleigh-Benard flow 

It is tempting to apply our considerations also to Rayleigh-Benard (RB) flow, which has 
recently been studied in detail [13-19]. RB flow is governed by a large scale "wind" 

^Strictly speaking, instead of the factor in brackets, one should write (1 + auik/ks)^^^"^)^^^^ , and a„ 
and as will not be independent. 
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T^ , |T6|, |T^, i.e., the temperature is subjected to a strong shear. This wind is of course 
caused by buoyancy and somehow a reminiscence of the convection rolls for lower Ra. Yet 
it may well be that the temperature in the interior can be considered as a passive scalar. 
This view is supported by the observation, that the local heat flux in the middle of the 



cell is tiny 16, 19 



First, we compare our predictions for the scaling exponents in shear flow with the 
experimentally measured ones and with those of the Bolgiano-Obukhov (B059) scaling 
theory , which has become the canonical way to account for the RB scaling exponents 
in the last years but which has also been criticized [^, |^, As we see from 



table 1, the experimental scaling exponents, with their limited accuracy, neither favor the 



shear flow nor the B059 explanation. Also the numerical simulations by Benzi et al. [|19 
though excluding a pure K41 scenario, may not be able to distinguish between the B059 
and the shear scenario. The numerical simulations by Kerr and Malevsky (3D) and 
by Werne (2D) do not favor B059. 

One can also speculate whether RB flow can be considered as helical flow (for a single 
realization due to spontaneous symmetry breaking) so that (|) can directly be applied. 
The effective forcing of the velocity field due to the heating will partly be helical so that, 
besides the energy fiux, a helicity fiux will be induced at large scales, resulting in a oc 
spectrum as discussed in Section 1. 

In the following we will analyze whether the crossover predictions of the two theories 
allow us to distinguish between them. Therefore we shortly repeat the main ideas of 
B059, which was originally suggested for stably stratified fiow [^. The kernel of the B059 
argument is dimensional analysis. Assuming that (3g and ee are the relevant parameters, 
where /3 and g are the volume expansion coefficient and the acceleration due to gravity, 
one obtains 0, gl|, 



for the velocity and temperature wave vector spectra. If, instead, ee and the kinetic energy 
dissipation rate e„ are the relevant parameters, one retrieves Kolmogorov's classical result 
(|l]). The crossover from B059 to K41 is expected for wave vectors k > ks, where the 
Bolgiano scale ks 0, ^ is obtained by equating (0) and (|^), 

kn = e^-'/W^\/3gf'. (7) 

With the exact result ^ ^ 

e„ = K^L-^PrRa{Nu - 1), e© = kA'^L-'^Nu, (8) 

the Bolgiano scale ks can be calculated. Here, Nu = H / {kAL~^) is the Nusselt number, 
H the heat fiux (= {u^Q) in the non diffusive case), Pr = u/k the Prandtl number, 
Ra = (3gAL^/ (uk) the Rayleigh number, L the height of the Rayleigh-Benard cell, A the 
temperature difference between the top and the bottom plate, z/ the viscosity, and k the 
diffusivity. For large Nu we obtain from eq. (0) k]^^ = Pr~^^^Ra~^^^Nu^^'^L. The K41 
range in the spectra should show up in between the B059 range and the viscous range, 
if k]j^ is larger than the typical crossover scale to the viscous subrange. For the velocity 
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structure function this scale is known |2^, ^ to be about lOrj, where rj = v^' e 



,3/4. -1/4 



Pr^/^Ra-^/^Nu-^I^L is the Kolmo gorov length |^ . For temperature structure functions 
and high Pr one could argue that k]^^ should be compared with IOt^o, where r^e = Pr~^/'^r] 



is the thermal Kolmogorov length which characterizes the crossover between the 

Batchelor regime (predicted for high Pr in a passive scalar intensity spectrum p9| ) 



and the VSR. For very low Prandtl number Pr~^/^ri should be used |29|. For existing 
experiments the resulting differences in the Pr dependence cannot be detected anyhow as 
Pr is always roughly 1, so this difference is only of theoretical interest for the moment. In 
what follows we will compare fc^^ with IOt] and give the Pr dependences resulting from 
a comparison of k]^^ with 10?7e (i.e., the large Pr case) in brackets. 

The condition A;^^ > 10?7 for the appearance of the K41 range in between the B059 
and the viscous range requires that 

Nu > ^ 22Pr (9) 

[22Pri/3 when using rj0 instead of r]]. I.e., beyond a certain crossover Rayleigh number 
Racoi defined by k]^^ = lOr/, there will be three ranges. For the Helium (Pr = 0.7) 
convection experiment [|1^, Nu = 0.165Pa^/^, or, embodying the suggested ^ Pr 
dependence, Nu = O.lSrPr-^/^Pa^/^. Thus, if the B059 scenario were true, pure B059 
scaling should be seen only for Ra < Raco ~ 3 ■ lO^Pr^ ^ 10^ [3 ■ lO^Pr^/^ r]e instead 
of 77], whereas for Ra > Raco ~ 10^ also the K41 range should be visible. Instead, pure 
Pe((x') oc ci;-i-3±o.i gcaiing jg observed at least up to Ra ~ 10^^ |T^. The discrepancy 



between Raco and 10^^ can only be resolved by allowing for some prefactors in the 
dimensional relations like eq. (|^) which are far beyond the order of 1. Note that Raco is 



near the experimentally observed transition between soft and hard turbulence |jT3|, |T^, . 
Our results predict that only for larger Pr the B059 scaling range might become more 
extended. Note that in this kind of analysis we neglected the height dependence of e and 
ee and thus of fcs by employing eqs. (|D0. 

For comparison, we now try to estimate the crossover between SDR and K41 in RB 
flow. This amounts to a comparison between the typical shear scale 

k7'=eu'/ys'/' (10) 

and lOrj. 

First, we estimate the shear s in RB flow. We do so by using the experimental infor- 
mation on the large scale velocity U = 0.16fi;L~^Pa^/^ [15| and on the mixing layer length 



Am, namely, s = U/Xm- Originally the mixing layer has been theoretically introduced [|I4| 
and later visualized ||30| as a region between the thermal boundary layer and the bulk of 

■^The Shraiman-Siggia theory ||2^, assumes Pr ^ 1, but it has proven to work also for smaUer Pr. 
Yet ah predictions on Pr dependences should be taken with care. 

^The nonuniversal spectral behavior for large oj and Ra > 10^^ jl^ could consistently be accounted 
for by considering the reduced high frequency resolution of the temperature probe because of its finite 
size [||. 

^In the numerical simulations of Benzi et al. a height dependence of ks is revealed. These 

simulations employ periodic boundary conditions. B059 scaling is found near the boundary (where 
is large), but in contrast to experiment not in the center of the cell. 
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turbulence. \m/L oc Ra'^f^ is found |T^, |T^. Embodying the predicted Pr dependence 
2^, this can be extended to \m/L oc Pr^/'^ Ra'^/'^ . In recent experiments []TB|, 
was redefined as distance of the maximum of the large scale velocity to the wall. In the 
employed L = 15.2cm RB cell |]T^, |T7|, a reduced high frequency resolution of the tem- 
perature probe (which is employed to measure A^, for the tricky experimental details, 
see ref. [^) should be expected around Rua = (lOOL^/a^)^/^ [^Sf, where a = 200/im is 
the extension of the probe. Indeed, for Ra < Raa the measured mixing length A^ 



agrees pTI with the suggested Ra scaling law Am ~ 1.2LPr^^'^ Ra ^1'^ . We can thus finally 



estimate 

s = U/X^ ^ OAKL'^Pr-^/'Ra'^/'^ (11) 

We now plug ([TT| ) and (||) into ( |T0|) and obtain ~ SLPr^^'^ Ra~^^'^^ . When com- 
paring to the crossover lOrj [or lOr^e] to the dissipation dominated range, we get the 
remarkable result, that the Ra scaling of and lOr] is the same, namely 

kj^/lOr] ^ O.SPr^/^ (12) 

[0.5Pr^/^ for rje instead of rj] for all Ra. I.e., for such small Pr as in the He cell (Pr = 0.7) 
the shear dominated range is directly followed by the dissipation range with no K41 range 
in between. For larger Pr the K41 range might be observable. 

Although the shear theory crossover scenario seems to agree better with experiment 
than the one for B059 (for small Pr = 0.7), many questions remain open and our analysis 
is far from conclusive, as it does not go beyond dimensional analysis. Moreover, the 
estimation s = f//Am in eq. (|ll]) can be criticized, because the shear might be less in 
the center of the cell where Pe(a;) is measured. Also, the z-dependences of the energy 
dissipation rate and of the shear s have been neglected. The reason for this section is 
only to demonstrate that the B059 scenario is not completely conclusive, either and that 
the shear scenario might well be able to account for the experimental observations. 



4 Summary 



To summarize, we derive the low frequency Pe{uj) oc uj~^^^ power law for passive scalars 
in a shear fiow [l^ . It should be followed by a K41 range (for high degree of turbulence) 
and of course by a viscous range. This scenario might also account for the measured 



temperature power spectra in RB fiow |jT5[, possibly without the K41 range in between 
SDR and VSR for low Pr. One of many open questions is the manner in which higher- 



order moments scale. Benzi et al.'s [T^ numerics show strong intermittent behavior, so 
dimensional analysis would fail here. We hope that this work motivates the analysis of 
shear effects for both the velocity and the temperature field in numerical simulations and 
in experiment. 
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Table 





K41 


B059 


shear flow 


experiment 




-5/3 
-5/3 


-11/5 = -2.2 
-7/5 = -1.4 


-7/3 = -2.33 
-4/3 = -1.33 


-2.2 ±0.2 p2| 
-1.3 ±0.1 [|15] 



Comparison of the scaling exponents predicted by different theories. The last column 
shows the experimental results. Tong and Shen's refer to velocity structure functions 
at Ra = 10^ — 10^°. I estimated the experimental error to be at least 0.2. Note that 
the very precisely measured Pe(ci;) scaling exponent [jl5| better agrees with the shear flow 
scenario than with B059. 



Figure 



Figure 1: Temperature power spectra in shear dominated flow, measured behind a wake 
of a heated cylinder. The data are taken from Sreenivasan's work [jlO[ with kind permission 
of the author. The spectrum shows a —4/3 power law. The dashed line shows the shear 
crossover frequency s. In the inset the power spectrum, multiplied by u^^^, is shown. 
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